Field-Induced Magnetic Ordering in the Quantum Spin System KCUCI3 
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KCUCI3 is a three-dimensional coupled spin-dimer sys- 
tem and has a singlet ground state with an excitation gap 
A/A;b = 31 K. High-field magnetization measurements for 
KCuCl3 have been performed in static magnetic fields of up 
to 30 T and in pulsed magnetic fields of up to 60 T. The 
entire magnetization curve including the saturation region 
was obtained at T = 1.3 K. From the analysis of the mag- 
netization curve, it was found that the exchange parameters 
determined from the dispersion relations of the magnetic ex- 
citations should be reduced, which suggests the importance 
of the renormalization effect in the magnetic excitations. The 
field-induced magnetic ordering accompanied by the cusplike 
minimum of the magnetization was observed as in the iso- 
morphous compound TICUCI3. The phase boundary was al- 
most independent of the field direction, and is represented by 
the power law. These results are consistent with the magnon 
Bose-Einstein condensation picture for field-induced magnetic 
ordering. 

PACS number 75.10.Jm 



I. INTRODUCTION 

Recently, the physics of coupled antiferromagnetic 
spin-dimer systems has been attracting considerable at- 
tention from the viewpoints of magnetic excitations 
p[-pjj| , the magnetization processes including magnetiza- 
tion plateaus |Tl|-|l8[] and field-induced three-dimensional 
(3D) magnetic ordering [^9|-p5[. This study is concerned 
with the magnetization process of and the field-induced 
3D magnetic ordering in KCUCI3, which is an S = h 
coupled antiferromagnetic spin-dimer system and has a 
singlet ground state with an excitation gap (spin gap) 

A/k B =31 K glUH. 

KCUCI3 has a monoclinic structure (space group 
The crystal structure is composed of pla- 



scattering [g,|,|30|-gj and ESR measurements |]2§. The 
effective dimer approximation |3],|34| and the cluster series 
expansion Q were successfully applied to analyze the dis- 
persion relations obtained. From these analyses, it was 
found that the origin of the spin gap is the strong anti- 
ferromagnetic interaction in the chemical dimer Cu 2 Clg, 
and that the neighboring dimers are weakly coupled not 
only along the double chain, but also in the (1,0,-2) 
plane, in which the hole orbitals of Cu 2+ spread. Con- 
sequently, KCUCI3 was characterized as a weakly and 
three-dimensionally coupled spin-dimer system. 

KCUCI3 differs from isostructural TICUCI3 and 
NH4CUCI3 in magnetic character. TlCuCh; also has 
a gapped ground state, however, the gap is consider- 
ably suppressed due to strong interdimer interactions 
§, ||, ||, 13 H]. On the other hand, NH 4 CuCl 3 has a 
gapless ground state at zero magnetic field |3(|. The 
remarkable feature of NH4CUCI3 is the existence of mag- 
netization plateaus at one-quarter and three-quarters of 
the saturation magnetization jl| ^7) . 

When a magnetic field is applied in the spin gap sys- 
tem, the gap A is suppressed and closes at H g = A/g^B- 
For H > H g the system can undergo magnetic order- 
ing due to 3D interactions with decreasing temperature. 
Such field-induced magnetic ordering was observed in 
the isostructural TlCuCl 3 |f-|5|]. For the field-induced 
magnetic ordering, two characteristic features which can- 
not be exp lained by the mean-field approach from the real 
space |38 ( [39) have been observed, irrespective of the ap- 
plied field direction, when the applied field H is slightly 
higher than H g . One is that the magnetization has a cus- 
plike minimum at the transition temperature Tn- The 
other is that the phase boundary between the paramag- 
netic phase and the ordered phase can be expressed by 
the power law 



P2 1 /c) )29|. 

nar dimers of CU2CI6. The dimers are stacked on top of 
one another to form infinite double chains parallel to the 
crystallographic a-axis. These double chains are located 
at the corners and center of the unit cell in the b — c plane 
as shown in Fig. 1. From the structural point of view, 
KCUCI3 was first assumed to be a double chain spin sys- 
tem magnetic excitations in KCUCI3 have 
been extensively investigated through neutron inelastic 
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with (j> = 2.1(1). Here, H^(T) denotes the transition 
field at temperature T. These features can be under- 
stood in terms of the Bose-Einstein condensation (BEC) 
of excited triplets (magnons) flll^-p^ . 

When H is slightly higher than H g , the density of 
magnon n which corresponds to the uniform magneti- 
zation is small. In this case, the system can be mapped 
onto the dilute magnon model. Using the Hartree-Fock 
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approximation, Nikuni et al. [|20| demonstrated that the 
magnetization has the cusplike minimum at the transi- 
tion temperature where the BEC of magnons occurs and 
that the phase boundary can be described by eq. (1) with 
4> = |. Although the theoretical exponent is somewhat 
smaller than the experimental value <f> — 2.1, the magnon 
BEC theory gives a good description of the experimental 
results for TlCuCl 3 . 

The magnetic anisotropy is negligible in KCUCI3 |27]] 
as it is in TICUCI3 [^,^4). Therefore, the previously 
mentioned characteristic features for the magnon BEC 
should also be observed in KCUCI3, because the magnon 
BEC is universal in the isotropic spin gap system. With 
this motivation, we carried out magnetization measure- 
ments in static high magnetic fields above the gap field 
H g 20 T using a hybrid magnet. Since the entire 
magnetization curve in KCUCI3 has not been observed, 
we also performed magnetization measurements in pulsed 
high magnetic fields of up to 60 T. As shown below, the 
saturation of the magnetization was observed at H s ~ 50 
T. From the saturation field H s , it was found that the ex- 
change parameters determined from the analyses of the 
dispersion should be reduced. This suggests the impor- 
tance of the renormalization effect in the magnetic exci- 
tations in KCuChj. 

II. EXPERIMENTAL DETAILS 

The single crystals of KCUCI3 were prepared by the 
Bridgman method. The details of sample preparation 
were almost the same as those for TICUCI3 p3fl . For 
KCUCI3, the temperature of the center of the furnace was 
set at 500°C. Single crystals of size 0.5~5 cm 3 were ob- 
tained. The crystals are easily cleaved along the (1, 0, 2) 
plane. The second cleavage plane is (0, 1,0). These cleav- 
age planes are perpendicular to each other. 

The high-field magnetization process in KCUCI3 was 
measured at T = 1.3 K using an induction method with a 
multilayer pulse magnet at the Research Center for Ma- 
terials Science at Extreme Conditions (KYOKUGEN), 
Osaka University. The magnetic fields were applied per- 
pendicularly to the (1,0,2) plane and along the [2,0,1] 
direction which is parallel to both cleavage (1,0,2) and 
(0, 1, 0) planes. We used a single crystal in the measure- 
ment for if||[2, Q, 1]. For _ff_L(l,0,2) measurement, we 
stacked several single crystals along the cylindrical axis 
of the sample holder. 

The temperature dependence of the magnetization of 
KCUCI3 was measured in static magnetic fields of up to 
30 T by a sample-extraction method with a 40 T-class 
hybrid magnet at the National Institute for Materials 
Science, Tsukuba. The magnetic fields were applied per- 
pendicularly to the cleavage (0, 1,0) and (1,0,2) planes, 
and along the [2, 0, 1] direction. These field directions 
were perpendicular to one another. 



III. RESULTS AND DISCUSSIONS 
A. Magnetization Process 

Figure 2 shows the magnetization curves of KCUCI3 
measured at T = 1.3 K for H _L (1, 0, 2) and H || [2, 0, 1]. 
Here, the magnetization curves are normalized with the 
fir-factors, i.e., g = 2.26 for H _L (1,0,2) and g = 2.04 
for H || [2,0, 1], which were obtained by ESR measure- 
ments. The magnetization curves for two different field 
directions almost coincide with each other when normal- 
ized by the ^-factor, although the data for H || [2,0,1] 
are somewhat scattered. This implies that the differ- 
ence between the transition fields H g and the satura- 
tion fields H s for H _L (1,0,2) and H || [2,0,1] is 
attributed to the anisotropy of the g-factor, and that 
the magnetic anisotropy is negligible, as previously con- 
cluded p7| ]. Due to the spin gap, the magnetization is 
almost zero up to the transition field (g/2)H g = 23 T, 
and increases rapidly and monotonically, and then satu- 
rates at (g/2)H s = 54 T. As shown in Fig. 2, the slope of 
the magnetization in the vicinity of H s and H s is steeper 
than that in the intermediate field region. This behavior 
arises due to the quantum fluctuation. 

In the magnetization slope region between H g and 
H s , transverse spin components have a long-range or- 
der. Since the lowest magnetic excitations occur at 
Q = (0, 0, 1) and its equivalent reciprocal lattice points, 
we assume that the spin ordering equivalent to that ob- 
served in TICUCI3 is realized in the ordered state (see 
Fig. 3). For the exchange interactions shown in Fig. 3, 
we use the notation given in ref . [Q : the main intradimer 
exchange is denoted as J. J(; mn ) and J'n mn \ denote the 
exchange interactions between dimers separated by a lat- 
tice vector la + mb+nc. Applying the mean- field approx- 
imation to the interdimer interactions, we calculate the 
magnetization curve for the ground state. The details of 
the calculation are given in the Appendix. 

In Table I, we list the exchange parameters in KCUCI3 
as determined by Mullcr and Mikeska H. They applied 
a cluster series expansion to analyze the dispersion re- 
lations observed in KCUCI3 ||. Their theory describes 
the experimental results, and the individual interdimer 
exchange parameters which cannot be obtained by the 
effective dimer approximation were determined. How- 
ever, they used the same intradimer exchange interac- 
tion, J — 4.25 meV, as that obtained from the effective 
dimer approximation (3). 

The dotted line in Fig. 4 is the magnetization curve 
calculated with the exchange parameters obtained by 
Miiller and Mikeska 0. We also plotted the experimen- 
tal result obtained at T = 1.3 K for H _L (1,0,2). The 
calculated gap field H g coincides with the experimental 
value, while the calculated saturation field H s is consid- 
erably larger than the experimental value. We infer that 
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the saturation field of eq. (A9) in the Appendix is close 
to the rigorous solution, because even for the S = | f D 
antiferromagnet, the classical calculation gives the same 
saturation field as the rigorous solution [[Io|. Therefore, 
we have to reduce the exchange parameters. However, it 
is difficult to adjust all of the exchange parameters, be- 
cause we have only two equations for H g and H s . Thus, 
we make a uniform reduction for all of the intcrdimcr 
exchange parameters. The thin solid line in Fig. 4 is 
the result calculated with the parameters listed in Ta- 
ble I. The intradimer interaction J and the interdimer 
interaction J(i mn ) or Jn mn \ are reduced by factors 0.9 
and 0.8, respectively. Both the calculated gap field H g 
and saturation field H s coincide with experimental val- 
ues. Hence, we suggest that the renormalization effect in 
the magnetic excitations cannot be neglected, although 
the present system is 3D. The renormalization factor n/ 2 
in the S = h ID antiferromagnet is well known [(41| . 
The difference between the calculated and experimental 
results for H g < H < H s is due to the quantum fluc- 
tuation, which is not taken into account in the present 
calculation. 

B. Field-induced magnetic ordering 

As mentioned in the Introduction, the present sys- 
tem can undergo magnetic ordering at a low tempera- 
ture, when the magnetic field H is higher than the gap 
field H g . Figure 5 shows the low-temperature magneti- 
zation measured at various external fields above H g for 
H || b, H _L (1,0,2) and H || [2,0,1]. With decreasing 
temperature, the magnetization exhibits a cusplike min- 
imum, irrespective of the field direction, as observed in 
isostructural TICUCI3. We assign the temperature with 
a cusplike minimum in the magnetization to the transi- 
tion temperature Tn. However, there is a certain degree 
of error in determining the transition points, because the 
experimental data are somewhat scattered. 

When the magnetic field is slightly higher than the 
gap field H g , the number of created triplets (magnons) is 
small. Through the transverse components of the inter- 
dimer exchange interactions, the magnons can hop to the 
neighboring dimers in the same way as particles, which 
have bosonic natures. The longitudinal component of the 
the interdimer exchange interaction gives rise to the in- 
teraction between magnons. Hence, the system can be 
mapped onto the interacting dilute boson system, and 
the phase transition can be described by the formation of 
the coherent state of dilute magnons, i.e., Bose-Einstein 
condensation (BEC) Jl9| pl|| . The cusplike minimum of 
the temperature variation of the magnetization is char- 
acteristic of the BEC of the dilute magnons |2^], and 
cannot be described by the mean-field approach from the 
real space |3^, It is known that the mean- field ap- 
proximation gives a reasonable description of the physical 



quantities at T = 0, while it does not at finite temper- 
atures. The mean-field approach gives the temperature- 
independent magnetization below Tn. The density of 
magnons n corresponds to the magnetization per Cu 2+ 
ion m, i.e., m — (g/2)fiBn. The reason that the magneti- 
zation increases again below Tn is because the increase of 
the number of condensed magnons is greater than the de- 
crease of the number of thermally excited noncondensed 
magnons, so that the total number of magnons increases 
below Tn- The present experimental result is consistent 
with the magnon BEC theory of the phase transition. 

The increase of the magnetization below Tn becomes 
smaller with increasing magnetic field. At H = 30 T, the 
density of magnons n is n ss 0.2. For this large value of n, 
the condition of dilute magnons is no longer satisfied. In 
the dense magnon region, the hopping of magnon is sig- 
nificantly suppressed, and thus, the mean-field approach 
from the real space may give a better description of the 
phase transition. 

The transition temperature Tn and transition field ob- 
tained for H || b, H _L (1, 0, 2) and H || [2, 0, 1] are plotted 
in Fig. 6. Since the phase boundary depends on the g- 
factor, we normalize the phase diagram by the g-factor. 
The g-factors used are g — 2.05 for H \\ b, g = 2.26 for 
H _L (1,0,2) and g = 2.04 for H || [2,0,1], which were 
determined by ESR measurements. Figure 7 shows the 
phase diagram normalized by the ^-factor. The phase 
boundaries for three different field directions almost co- 
incide when normalized by the g-factor, although the 
boundary for H _L (1,0,2) tends to deviate from the oth- 
ers at around (g/2)H ~ 30 T. This result may reconfirm 
that the magnetic anisotropy is negligible in KCUCI3 , as 
in TlCuCl 3 . 

The phase boundary near the gap field can be de- 
scribed by the power law of eq. (1), as predicted by 
the magnon BEC theory pq] . We fit eq. (1) to the data 
for T < 3 K, for which n < 0.03. The solid line in Fig. 
7 is the fitting with 4> = 2.3(1) and (g/2)H s = 22.6(1) 
T. This exponent cf> is close to <f> = 2.1(1) observed in 
TlCuCl 3 ||,||]. The exponent (j> = 2.3 obtained by the 
present measurements is somewhat larger than the value 
0=1 predicted by the magnon BEC theory [po| . We 
note that the exponent depends on the temperature re- 
gion used for fitting. If we use the data up to higher 
temperatures, we have a larger exponent. This is be- 
cause the ordered phase exists in a closed area in the 
magnetic field vs temperature diagram, and thus the en- 
tire phase boundary cannot be described by the power 
law with a single exponent. Since the low-temperature 
data contribute greatly to the determination of the ex- 
ponent 4> for the phase boundary near the gap field, the 
measurements at lower temperatures are needed. 

The increase of the transition field obeying the power 
law with increasing temperature cannot be described in 
terms of the mean-field approach from the real space 
p8],|3!l or the temperature dependence of the excitation 
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gap |33|. The mean-field result gives the transition field 
almost independent of temperature at low temperatures. 
The temperature dependence of the excitation gap was 
investigated by means of neutron scattering ]3^ |. The 
gap is almost independent of temperature up to approx- 
imately 10 K, and then increases significantly. 

IV. CONCLUSION 

We have presented the results of magnetization mea- 
surements on the coupled spin-dimer system KCUCI3. 
The entire magnetization curve including the saturation 
region was obtained at T = 1.3 K, using pulsed magnetic 
fields up to 60 T. Applying mean-field approximation on 
the interdimer interactions, we analyzed the magnetiza- 
tion curve. It was found that the exchange parameters 
determined from the dispersion relations of the magnetic 
excitations should be reduced to fit both the gap field 
and the saturation field. This suggests the importance of 
the renormalization effect in the magnetic excitations. 

The field-induced magnetic ordering was investigated 
in static magnetic fields up to 30 T. As observed in 
isostructural TICUCI3, the magnetization exhibits a cus- 
plike minimum at the transition temperature. The phase 
boundary is almost independent of field direction, and is 
described by the power law. These features are compati- 
ble with the magnon BEC theory. Therefore, we conclude 
that the field-induced 3D magnetic ordering near the gap 
field in isotropic spin gap systems such as KCUCI3 and 
TlCuChj is universally represented by the BEC of dilute 
magnons. 



structure in the ordered phase in KCUCI3 is equivalent to 
that observed in TICUCI3 (see Fig. 3) p5[ . We express 
the spin state of a dimer with the total spin S = S\ + S2 
and the z component S z = Sf + S% as \S,S Z ), where 
subscript numbers 1 and 2 distinguish spins in a dimer. 
In the ordered phase, the transverse components of two 
spins on a dimer should be antiparallel and their magni- 
tudes should be the same, i.e., (S^' v ) = — (S%' y ). This 
indicates that the spin state |1, 0) does not contribute to 
the basis state. Thus, we write the basis state of the j-th 
dimer as 

ipj = |0,0) cos (9+ (|l,l)cosy> e 1 ^ - |l,-l)sin<p e _i< ^) sin (9 . 

(Al) 

Angles 9 and <p were introduced to satisfy the normaliza- 
tion condition. The phase (f>j corresponds to the angle 
between the x-direction and the transverse component of 
the spin on the j-th dimer. 

The average values of spin operators are given by 

(S],i) =(S|, 2 ) = isin 2 flcos2^, 

(Sf A ) = - (S+ 2 ) = — -j= cos 6 sin (cos (p + sin <p) e — , 

(S^) = - (Sr 2 ) = — -= cos 9 sin 9 (cos ip. + sin ip) . (A2) 
V ^ 

From the spin structure shown in Fig. 3, we see that 
all of the phases <f>j of dimers are the same in a chemical 
double chain, and that the phases of dimers in the double 
chains located at the corner and the center of the unit 
cell in the b — c plane differ by ir. With eq. (A2) and the 
exchange interactions shown in Fig. 3, the energy per 
dimer is expressed as 
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APPENDIX A: DERIVATION OF THE 
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MAGNETIZATION CURVE cff • ; 



i - if" J 1 " 111111 

We calculate the magnetization curve at T = 0, as- reff — { t V 

suming that the triplet excitations are created only on ^ ,! ' 2 ^ I \ \ 'i'*) l>; 

the dimer site. The interdimer interactions are treated T cS _ T i t 

by the mean-field approximation. In the magnetization ' ' 2 ( ' ' 

slope region, the transverse spin components form the an( j 
long-range order. Since the lowest magnetic excitations 

in KCUCI3 and TICUCI3 occur at the same reciprocal lat- 7 — T -l — T' -l t j_ j 1 j- — t' 

tice points, i.e., Q = (h,0J) with integer h and odd / in J ~ J{W0) + 2 J ( W ^ +J (^U) + J (i,|,|) + 2 J ( 2 ' ' 1 ) ' 

the a* — c* plane [0,|l||,pf , we can assume that the spin (A6) 



4 



Angles 9 and <p can be determined by dE/dO = and 
dE/dip = 0, which lead to 

J — gfJ-sH cos2</3 + J cos 26 (1 + sin 2^) + J sin 2 9 cos 2 2y> = 

(A7) 

and 

g/iB-ff sin 2ip + J cos 2 cos 2ip ~ J sin 2 cos 2y> sin 2<p = , 

(A8) 

for the ordered state between the gap field H g and the 
saturation field H s . The magnetization curve is obtained 
by solving eqs. (A7) and (A8) self-consistently. 

At the saturation field H s , sin 6 = 1 and simp = 0. 
Thus, H s is given by 

gfi B H s = J — J + J . (A9) 

At the gap field H g , s'md = 0. Substituting this condition 
into eqs. (A7) and (A8), we obtain 

g^ B H g = V J 2 + 2J.J . (AiO) 

It is noted that the right-hand side of eq. (A10) is equiva- 
lent to the lowest excitation energy given by the effective 
dimer approximation, in which the individual interdimcr 
interactions J(i m n) and J'(i mn ^ are reduced to an effective 
interdimer interaction J^ nn \ as shown in eq. (A5) ||. If 
we neglect the highest |1,— I) state, i.e., sinip = 0, we 
obtain 

gf x B H s = J + j. (All) 

This result can also be derived from the Tachiki-Yamada 
theory |pq,|39||. 
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TABLE I. Exchange interactions in KCUCI3 in the unit of 
meV. 



J 


^(lOO) 


J (100) 


J (i>U) 




"^(2,0,1) 


ref. 


4.25 
3.83 


-0.021 
-0.017 


0.425 
0.340 


0.850 
0.680 


0.170 
0.136 


0.799 
0.639 


this work 



G 



FIG. 1. Crystal structure of KCUCI3 viewed along the 
a-axis. Shaded, small open and large open circles denote K + , 
Cu 2+ and CI" ions, respectively. 
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FIG. 2. Magnetization curve of KCUCI3 measured at 
T — 1.3 K for the magnetic fields H || [2,0,1] and 
H _L (1, 0, 2). The values of the magnetization and the mag- 
netic field are normalized by the gr-factor. 



FIG. 3. Important exchange interactions in KCuCl3 and 
the spin structure observed in the ordered phase of TICUCI3 
for H || b. The double chain located at the corner and the 
center of the chemical unit cell in the b — c plane are repre- 
sented by solid and dashed lines, respectively. The shaded 
area is the chemical unit cell in the a — c plane. 
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FIG. 4. Magnetization curve calculated by the mean-field 
approximation. The dotted and thin solid lines are the 
results with the exchange parameters obtained by Miiller 
and Mikeska and reduced ones, respectively. The thick 
solid line is the magnetization curve observed in KCUCI3 for 
HI. (1,0,2). 
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(b) 



FIG. 5. Low-temperature magnetizations of KCUCI3 mea- 
sured at various external fields for (a) H \\ b, (b) H ± (1, 0, 2) 
and (c) H || [2,0,1]. 
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FIG. 6. Phase boundaries in KCUCI3 obtained for three 
different field directions. 
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FIG. 7. Phase diagram in KCUCI3 normalized by the 
p-factor. The solid line denotes the fitting by eq. (1) with 



(g/2)H s = 22.6(1) T and 
visual guide. 
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